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Abstract 

The thermodynamic Bethe ansatz is applied to a quantum integrable 
spin chain associated with the Lie superalgebra osp{l\2). Using the string 
hypothesis, we derive a set of infinite number of non-linear integral equa- 
tions (thermodynamic Bethe ansatz equation) , which characterize the free 
energy. The low temperature limit of the free energy is also discussed. 
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Solvable lattice models related to Lie superalgebras [Q have received much at- 
tentions |[ ^ ||, ||, 0, ^, ^. Characteristically, these models possess both 
fermionic and bosonic degree of freedom and their Boltzmann weights {R- 
matrices) satisfy the graded Yang-Baxter equations Q. Some of them include 
well-known models as special cases. For example, supersymmetric t — J model 
is related to the Lie superalgebra s^(l|2). To analyze such models, Bethe ansatz 
is often used (See for example [| |^, |l2|, g ^ |l^, ^ |l9|, |2^, 



2^ |2^, ^ |2^, |2^, ^ and references therein). Moreover, thermodynamics 
have been discussed by several people. In particular, the thermodynamic Bethe 
ansatz (TBA) equation was given for the supersymmetric t — J model [|o[ ^ 
and the supersymmetric extended Hubbard model which is related to 

sl{2\2). More recently, the TBA equation for sl{r\s) model was discussed 
in relation with the continuum limit of integrable spin chains. 

However, in contrast to sl{r\s) case, the thermodynamics of quantum spin 
chain related to the orthosymplectic Lie superalgebra osp{r\2s) is not much un- 
derstood. In particular, to the author's knowledge, there is no literature on the 
thermodynamic Bethe ansatz (TBA) equation even for the simplest orthosym- 
plectic osp(l|2) integrable spin chain |, |, |, 0, |, |, ^, [l7| . This is regrettable 
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because this model may be related to interesting problems such as = 1 super- 
symmetry in particle physics, and the loop model, which will describe statistical 
properties of polymers in condensed matter physics [ p7[ . 

The purpose of this letter is to construct the TBA equation for osp{l\2) 
integrable spin chain based on the string hypothesis |3^]. The resultant 

TBA equation, which describes the free energy is a set of infinite number of 
non-linear integral equations. We have also confirmed the fact that our TBA 
equation coincides with the one from the quantum transfer matrix (QTM) 



approach 



35| based on the transfer matrix functional relations [psj , 
nj, |23|, |2j (osp(l|2) version of the T-system). This letter yields a basis 
of future studies of thermodynamics for more general models such as osp{r\2s) 
model. 

Now we shall describe the ospil\2) model|, |, |, 0, |, |, |l6[ [l^. The R- 
matrix of this model is given as a rational solution of the graded Yang-Baxter 
equation associated with the three dimensional representation of the Lie 
superalgebra osp{l\2) (to be precise, 'super Yangian' Y{osp{l\2))), whose basis 
is Z2-graded and labeled by the parameters p{l) = p{'3) = l,p{2) = 0. One can 
construct the i?(u)-matrix as a 'Baxterization' of the Temperley-Lieb generator 
It reads explicitly as follows: 

R{u)=I + uPa (1) 

where denotes the graded permutation operator which has 9x9 matrix ele- 
ments {P^)ai = (— l)^^°^^^''^^a,d^6,c; -E^ is a 9 X 9 matrix whose matrix elements 
are [E^Y^ = aabO:~J; a is a 3 x 3 matrix: 




a = 1 . (2) 



There are two kinds of the definition of the transfer matrix for each model related 
to Lie superalgebras. The one is defined as the supertrace of a monodromy 
matrix and the other is the ordinary trace of a monodromy matrix. In this letter, 
we adopt the latter. In this case, the transfer matrix T{u) of the corresponding 
system can be written as 

T{u) = tra [RaN [iu) ■ ■ ■ Rai (iu)] , (3) 

where N is the number of lattice sites; a denotes the auxiliary space; Raj{iu) 
denotes R(iu) which acts non-trivially on the auxiliary space and j-th site of the 
quantum space; R{u) — PR{u): Pf^^ = 6a,dSb.c is the (non-graded) permutation 
operator. 

The Hamiltonian is defined by taking the logarithmic derivative of above 
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transfer matrix at u = 0, 



i du 



u=0 



= ^Efe+i + ^^.+ij. (4) 



where we assume a periodic boundary condition. Here J is a real coupling 
constant which characterizes the phase of the model; the ferromagnetic and 
antiferromagnetic regimes correspond J > and J < 0, respectively. The 
eigenvalue formula of (||) (also written as T{u)) is given ^ |l^, ^ in the dressed 
vacuum form as follows 

Tin) = i-ir~-iu+r^^!^+u- 




Q{u+^)Q{u + i) 

(■ \ ^ 
\ Qiu + 2i) 

where Q{u) = rij=i(" ~ ""i); ^ {0, 1, ■ • • , N} is a quantum number; Uj G C. 
The dressed vacuum form is built on the pseudo-vacuum state corresponding 
n = 0. This formula (H) is free of poles under the following Bethe ansatz 
equation (BAE) 

= -(-1)--" ^^:!!^;^^+^! for ..{1,2,...,.}. (6) 
Q{uk + ^)Q(uk - i) 

Taking the logarithmic derivative of (^), we obtain the eigenvalue formula of 
the Hamiltonian (^) 

Following Takahashi's argument (see also ||36|]), we shall discuss the ther- 
modynamics based on the string hypothesis. The m-string solution of the BAE 
(H) is given |l^ as follows 

<'"-4" + ^(™ + l-2a), (8) 

where are the center of the strings; a G {1,2,..., m}; fc S {1,2,..., «,„}; 
rim is the number of the m-strings. Assuming (^) and taking the product of (||) 
on each string, we obtain 




TO 



N 



= (-i)(^-"+'^'"nn^™'«-4). (9) 

1=1 3 = 1 
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where e{u) := (u + ^)/{u — |) and 



Emiiu) = e{- :]e'{- ^ •••e 



\l-m\J V|;-m| + 2/ \l + m-2j \l + m 

-11 \ -1 f ^ \ -1 ^ ^ 



xe" ^ , T e-M . r •••e-M- if m ^ /, 



Taking the logarithm of we obtain 



^U27r/r + 5:^e,„,(ur-«^), (11) 
' 1 1=1 



1=1 j=l 

where 6{u) = 2arctan(2u), 7™ G and 
Qmiiu) = o(-^)+2e(- -]+... + 29 



\l-m\J \\l-m\ + 2j \l + m-2 



if 771 ^ I, 



2ei^\+2e[^ + ... + 2e' " 



21 V4/ \27n-2J V2m/ 

Substituting (||) into (Q), we obtain the energy in terms of the centers of the 
strings: 

{oo nm I 
EE 5^--)- («) 

Now we shaU take the thermodynamic hmit {N ^ oo) of (|ll|) and (^3|). Let 
particle, hole and vacancy densities be pf„(u), p^yi{u) and pm{u), respectively: 

jm jrn 
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In the thermodynamic hmit, ( |ll| ) reduces to the fohowing integral equations: 

oo 

!m{u)=pt{u)+Y^Almp'i{u), (15) 



1 = 1 

where ^/,„ is defined as follows 



Al„^ = An; = [|^ - m\] + 2[\l - m\ + 2] + ■ ■ ■ + 2[l + m - 2] 
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-m\ + l]-[\l-m\+3] [l + m-1], (16) 

and [m] acts on any function g{u) as follows 

[m]g{u) = /„i * g{u) for m 7^ 0, 
[Q\g{u) = 5(w), (17) 

where * denotes a convolution 

/oo 
f,niu- v)g{v)dv, (18) 

and fm{u) is defined as follows 

777- 

The operators ( [T^ ) and ( p7|) a rc osp(l|2) version of the ones introduced by 
Takahashi M. The energy (|13|) per site becomes 



iV 



The entropy per site is given as follows 

JJ - ^bY. {{pt{u) + pl,{u))Hpt[u) + p^^{u)) 

m=l 

-P:'„ (u) In p:^. (^.) - (u) In pP, (z.) } du, (21) 



where fcs is the Boltzmann constant. In the thermodynamic limit, the free 
energy F = E — TS (T: temperature) should be minimized with respect to 
density functions (|l^. Namely, we have 



m=l " V ^ ' 
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One can derive the following relation from (p^. 

OQ 

-^p^M^-E^'-^^/^rH- (23) 

1=1 

Combining (|2^) with (p^, we obtain the following integral equation 

oo 

ln(l + 7?™(m)) = 2tt/3JU{u) + E^™' 1^(1 + Viiu)~^), (24) 



1=1 



where P = l/ksT and rym(w) := pJ'„(u)/pP„(u). 

To proceed further, we assume the function r]m{u) (m > 1) have constant 
asymptotics in the limits u — > oo. After performing the Fourier transformation 
on both sides of above equation (^) , we obtain 

oo 

.F[ln(l + r7„0](fc) - 2nJ(3e-'^^''^+J2^[Ami]{kMhi{l + rjr^)]{k), (25) 

1=1 

where fc e M and 



^ ' smh ^ 

In the above, we adopt 

Him = / f{u)e''''^du, (27) 

-/ — OO 

as the Fourier components of a function f{u). Multiplying both sides in ( |25| ) by 

^""^^^"'^"'"^ {C™+!+5n,™-i form>2, ^2^) 

and taking summation with respect to m, we obtain 

T[\nl^^][k) = -^__l-^(^[ln(l+,^-i)](fc)-.F[ln(l+r;2)](fc)), 

T[\Tirim]{k) - ^^(^[ln(l + r;„,_i)](fc)+.F[ln(l + r;„,+i)](fc) 
2 cosh I 

-J^[ln(l + ?7,7/)](fc)) for TO > 2. (29) 
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By performing the inverse Fourier transformation, we obtain the foUowing TBA 
equation 

Inryi(w) = X * ln(l + 7?r^)(u) + X * ln(l + 772)(w), 

cosh TTU 

ln?7„(u) = K*\n{l+T]„,_i){u)-K*ln{l + r]-^^){u) 

+X*hi(l + 77,„+i)(-it) for TO > 2, (30) 

where the integral kernel K(u) is defined by 

K(u) = — . (31) 

^ ^ 2cosh7rw ^ ^ 

After performing the Fourier transformation on (^|) and multiplying by _Bj^ , 
we get 

= ^^^2 a 1 {} - 2cosh|.FK](fc) +.F[p^](fc)j , 
np'jik) = - 2cosh!i-i (2cosh^^[p:jj(fc) - npt-iKk) np'UiKk?! 

for TO > 2. (32) 

Combining above equations with (|2^) and (|l]) , and using the fact that r]„i (u) = 
Vm{—u) for TO > 1, we can derive the free energy per site / as 

f ^ J (^-^-l^ -keT j R{u)\n{l + r]i{u))du, (33) 

where R{u) is written as 

2sinh^ ^^^^ 
\/3 sinh 27rit 

To determine the unique solutions of the TBA equation (^) , we need to obtain 
the asymptotic values 77m(±oo) :— lim„^±oo f7m('") {m > 1)- To determine 
r7m(ioo), we shall consider the high temperature limit T — s- oo. In this limit, we 
assume that all the functions rjmiu) (to > 1) are independent of the parameter 
u. Then (30) can be rewritten as follows. 

2 1 + ^/2 

Vi - 



= + for m>2. (35) 

1 + Vrn 



The free energy ( |33| ) must have the following form. 



lim = -ln3, (36) 
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which follows from the fact that the freedom of the state per site is three, namely 
the entropy per site must be 5 = fcs In 3 in the high temperature limit. Thus 
we have rji — 2 and the following solutions of (p5|): 



m(m + 3) ^ _ , , 

?7m = -^—^ for TO > 1. (37) 

Therefore at any finite temperature (T > 0), the asymptotic values of r/rn(u) 
( |30| ) have the same values as (|37|): 

lim ry^^I^^i!^ forTO>l, (38) 



since 77„j(±cx3) also satisfy the same equations as (35) and assumed to be inde 



pendent of T . Consequently, the free energy per site at any finite temperatures 



can be uniquely determined by (p3|) through the TBA equation (30) and the 



asymptotic values (38). The above representation for the TBA equation (| 



the asymptotic values (38) and the free energy ( |33| ) can also be confirmed by 
the QTM method, which is free from the string hypothesis One might sus- 
pect that our TBA equation ( ^o|) is nothing but a certain specialization of the 
one js^ for the Izergin-Korepin model since the affine version of osp(l|2), i.e. 
osp(l|2)'^^' has close resemblance to the affine Lie algebra J^-f^ ■ However, to the 
author's knowledge, so far the TBA equation for the Izergin-Korepin model has 
been constructed |3^ so that it reduces to the one for a sM(3)-invariant model 
in the isotropic limit. Therefore, we expect that our TBA equation ( ^0| ) is new. 

From the above expression, one can analyze the the low temperature limit 
of the free energy. Let £,„(") = kBT\nri„i{u). In the case J > 0, we can see 
that the functions Em (to > 1) are always positive by applying the iteration 



method to the TBA equation (|3C|). Therefore from (24), the functions em{u) 
are written, in the limit T ^ 0, as 

e™(u) = 27rJ/„i(w) for m > 1. (39) 

Thus we have 

T],n = OO, P?n = 0. (40) 

These results mean that all the spins are in the same direction, which agrees that 
the J > case corresponds to the ferromagnetic regime. From above Eq. ( pO[ ) 
and Eq. (|20|), we have the ground state energy per site cq: 

eo = - J. (41) 

While in the case J < 0, one sees from the TBA equation ( |30|) that ry„i < 
to(to-|-3)/2 (m > 1). Therefore the second term in the equation (p3) will vanish 
in the limit T ^ 0. Thus the ground state energy per site cq is given by 

/ 47r 



eo = J - 1 ) , (42) 



which coincides with the result in fVi 
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